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Abstract 

In this paper we use a constructive approach based on gauge invariant description 
of massive high spin particles for investigation of possible interactions of massive spin 
2 particle. We work with general case of massive spin 2 particle living in constant 
curvature (A)dSd background, which allows us carefully consider all flat space, massless 
or partially massless limits. In the linear approximation (cubic terms with no more than 
two derivatives in the Lagrangians and linear terms with no more than one derivative in 
gauge transformations) we investigate possible self-interaction, interaction with matter 
(i.e. spin 0, 1 and 1/2 particles) and interaction with gravity. 
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Introduction 



In all investigations of massless particles interactions gauge invariance plays a crucial role. 
Not only it determines a kinematic structure of free theory and guarantees a right number 
of physical degrees of freedom, but also to a large extent fixes possible interactions of such 
particles. This leads, in particular, to formulation of so called constructive approach for 
investigation of massless particles models [TJ [21 El ffl E]- In such approach, starting with 
appropriate collection of free massless particles and requiring conservation of (modified) 
gauge invariance in the process of switching on an interaction, one can consistently reproduce 
such physically important theories as Yang-Mills, gravity and supergravity. 

Usual description of massive particles does not include gauge invariance and it is hard 
to formulate one simple principle for constructing consistent theories with such particles. 
A lot of different requirements, such as conservation of right number of physical degrees 
of freedom, smooth massless limit, tree level unitarity and causality, was used in the past 
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There exist two classes of consistent models for massive high spin particles, namely, for 
massive non-abelian spin 1 particles and for massive spin 3/2 ones. In both cases masses of 
gauge fields appear as a result of spontaneous gauge symmetry breaking. One of the main 
ingredients of this mechanism is the appearance of Goldstone particles with non-homogeneous 
gauge transformations. This, in turn, leads to the gauge invariant description of such massive 
spin 1 and spin 3/2 particles. But such gauge invariant description of massive particles could 
be constructed for higher spins as well, e.g. [Ill [13 [16l [13 [111 [III [20l [211 [221 ^ This allows 
one to extend the above mentioned constructive approach to the case of massive high spin 
particles. 

In the first section of our paper we give two simple examples of such constructive 
approach. One of them deals with triplet of massive spin 1 particles with gauge group 
SU(2) ~ 0(3). We show that constructive approach allows one to reproduce two well known 
possibilities: one based on the non- linear cx-model [23] , and, with the help of introduction of 
additional scalar field, usual model of spontaneous symmetry breaking with the doublet of 
Higgs fields. Another example devoted to electromagnetic interaction for massive spin 3/2 
particle j9j [10]. In the first linear approximation (i.e. cubic terms in the Lagrangian and 
linear terms in the gauge transformations) we construct the most general gauge invariant 
Lagrangian. Our results unambiguously show that any model with minimal gauge interac- 
tion of massive spin 3/2 particle in flat space must be a part of some (spontaneously broken) 
supergravity theory. 

Main part of our paper devoted to massive spin 2 particle interactions. It is an old but 
still unsolved problem and the main difficulty is connected with the so called sixth degree 
of freedom |25j . Indeed, for the free particle it is easy to choose a structure of mass terms 
so that there are exactly five degrees of freedom in the model. But in the interacting theory 
this dangerous (because it is a ghost) sixth degree of freedom reappear. Note that it is, in 
principle, possible that in the full non-linear theory this DoF become physical [2S112Z], but in 
what follows we will not consider such possibility. One more difficulty that appears in such 
theories is the absence of smooth massless limit [23 129]. Moreover, when one considers such 
theory in a background space with non-zero cosmological term (de Sitter or anti de Sitter 
spaces) one faces an ambiguity between flat space and massless limits [3"0~ 1 I3"TI 132] . 
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Here we use constructive approach based on the gauge invariant description of massive 
particles for the investigation of possible massive spin 2 interactions. It turns out that 
some results depends on the dimensionality of space-time as well as on the presence (or 
absence) of cosmological term. So we start in section 2 with the gauge invariant formulation 
of free massive spin 2 particle in the space-time of arbitrary d > 2 dimension with non- 
zero cosmological term, that could be positive or negative. This includes, in particular, the 
so called partially massless spin 2 particle, which could exist in the de Sitter space only 

[331 Ell sails Egg. 

Then, in section 3 we consider possible self-interaction of such particles in linear approx- 
imation. We will not make any suggestion about gauge algebra which could stand behind 
such theory nor we will not insist on any possible geometrical interpretation. Instead, we 
use a "brute force" method for our construction. Namely, we write the most general Lorentz 
invariant cubic vertexes in the Lagrangian as well as the most general ansatz for gauge trans- 
formations linear in fields and require that the Lagrangian will be gauge invariant. The only 
restriction (besides Lorentz invariance) we impose is the restriction on the number of deriva- 
tives. In this paper we consider only interactions with no more than two derivatives (and 
gauge transformations with no more than one derivative) leaving investigation of possible 
higher derivatives interactions for the future. 

Clearly, the existence of linear approximation does not guarantee that the construction of 
full non-linear theory is possible because obstruction can appear in the next approximations. 
Nevertheless, the investigation of linear approximation is a very important step, because 
structure of this approximation (and its whole existence) does not depend on the presence 
of any other field in the theory. Only in the next quadratic approximation one faces the 
problem that the closure of gauge algebra requires the presence of right (finite or infinite) 
collection of fields. Thus, the results obtained in our paper are essentially model independent 
(up to restriction on the number of derivatives). 

For the generic values of mass and cosmological constant the linear approximation for 
self- interaction exists in any space-time dimension d > 3, but for partially massless spin 
2 particle it exists in d = 4 dimensions only. This clearly related with the fact that only 
in d — 4 partially massless theories are conformally invariant [37]. One of the non-trivial 
results is that the structure of gauge transformations for Goldstone field turns out to be 
non-canonical, so if one tries to interpret (as usual in gravity) the £ M gauge transformations 
as general coordinate ones, then this field should give some non-linear realization of this 
transformations. One more non-trivial result is that two gauge transformations — vector 
£ M and scalar A ones, do not commute when the mass m^fl forming an algebra similar to 
Weyl one. 

In the next section we consider possible cubic couplings of massive spin two particles 
with matter. Namely, we consider interaction with (massive or massless) spin 0, spin 1 and 
spin 1/2 particles. In all cases (except massless spin 1 in d — 4) our results show an am- 
biguity between flat and massless limits, which shows itself through non-trivial dependence 
of coupling constant for scalar Goldstone field tp and matter fields interactions on mass and 
cosmological constant. As a result, partially massless spin 2 particle can interact with mat- 
ter having traceless energy-momentum tensor only, the most important example being the 
coupling of partially massless spin 2 with massless spin 1. 

Our last section deals with possible interactions for massive spin 2 particles with gravity. 
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As is well known, it is impossible to construct non-trivial consistent interacting theory for 
collection of massless spin 2 particles [38j ESI HQl EH H2], but it still leaves the possibility 
to construct non-trivial interaction of gravity with (one or many) massive spin 2 particles. 
Indeed, such possibility attend much interest recently in very different contexts. First of all, 
we show that it is impossible to construct cubic vertex with two massless and one massive 
spin 2 particles while vertex with two massive and one massless ones does exist. Moreover, it 
is easy to construct general covariant vertex with two massive fields and arbitrary number of 
massless gravitons. As for the interaction of partially massless spin 2 particles with gravity, 
no restrictions on the space-time dimension arise in this approximation. 



1 Toy models 

In this section we present two simple examples of constructive approach to massive high 
spin particles interactions. First one deals with the triplet of massive spin 1 particles, while 
second one devoted to possible gauge interactions for massive spin 3/2 particle. Results 
of this section where already known before, but it is instructive to see how they can be 
reproduced using gauge invariant formulation alone. 



1.1 Spin 1 

In this subsection we illustrate our approach by the simplest non-abelian 577(2) ~ 0(3) 
gauge theory with triplet of vector fields A^ a , a = 1,2, 3. It is well known that with the help 
of triplet scalar Goldstone fields ip a one can easily construct gauge invariant description of 
free massive particles. In this, the Lagrangian and gauge transformations look like: 

1 Tft^ 1 

4> = --{A^f + —{A^f-mA^d^ + -{d^f 

w = ^r, <w a = ™r (i) 

Now, if we switch on usual non-abelian gauge interaction: 

A, u a = d^AS - ge^AjAS - fa » v) 

5oV = d,C-ge abc A, b C (2) 

then the Lagrangian will not be gauge invariant any more: 

5 C = -gme abc A^d^ b C (3) 

but at the linear approximation the gauge invariance could be easily restored with additional 
vertex in the Lagrangian and appropriate corrections to gauge transformations: 

A = axe-fc^y V> s i<P° = Si£ a Vr (4) 

provided a\ = —g/2, g± = —g/2\. Note, that non-canonical value g/2'va. the gauge trans- 
formations of scalar fields ip a shows that these fields do not transform under linear triplet 
representation of gauge group. As we already mentioned in the introduction, this first step 
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does not depend on the presence of any other fields in the model. Now we can proceed in at 
least two different ways. At first, we can evade introduction of any additional fields and go 
on by introducing into the Lagrangian terms of the type A^ip 2 d^ip, (p 2 (d^(p) 2 and so on. It 
would lead us to the essentially non-linear theory with the Lagrangian 



C = — A {A, u a f + —(V) 2 - mA«E ak {ip)d^ k + ^WM^ (5) 

and gauge transformations for scalar fields: 

5<p k = m{E- 1 ) ka i a (6) 

where g kl = E ak E al . This Lagrangian will be gauge invariant provided the "triad" E ak 
satisfy the equation 

It is easy to check that non-canonical value g/2 in the linear approximation for the ip a gauge 
transformations is a direct consequence of this equation. 

But there is another way. We can introduce one more scalar field \ and try to stop 
iterations at some order. And indeed, if we add to the Lagrangian of linear approximation 
all possible cubic and quartic terms with new field: 

C 2 = \{d, X ) 2 + a 2 x(cML) V + %X W + a,m X {A, a f + a 5 ( V) V + a eX 2 ( V) 2 (8) 

as well as appropriate corrections to the gauge transformations: 

8X = 92(<P£), 5 2 ^ = g 3X C (9) 
we can obtain full gauge invariant theory, provided: 

9 _ _ 9 _ g 2 

2> a 3 — 9-, 93 — a 4 — — ^ j a 5 — a 6 — -g- 



2m 

iiiiiieu vtiiifciuit; X ~ X — 

rewritten in a familiar form: 



Moreover, if we introduce shifted variable x = X ~ then total Lagrangian could be 



C = -\{A^) 2 + \{d^Y + \{d,x) 2 -\e abc A^ 



+ | V(X^° - d,XV a ) + v( V) 2 (^ 2 + X 2 ) (10) 



Up to the arbitrary potential depending on the invariant combination if 2 + x 2 only this is 
just the well known model for spontaneous gauge symmetry breaking through the Higgs 
mechanism. In this, four scalar fields (p a and x are transformed under the linear doublet 
representation of gauge group. 

Note, that one more possibility arises when one consider systems with infinite number of 
gauge fields such as so called Kaluza-Klein models e.g. [121 HH SSI SSI 113 HB]- Indeed, let 
us consider simplest example — five-dimensional Yang- Mills theory with gauge fields A^, 
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where M = 0,1,2,3,5, and some non-abelian gauge algebra [t a ,t b ] = f abc t c . Suppose now 
that one dimension is a compact one, being a circle of radius R. From a four- dimensional 
point of view, such fields represent combinations of vector and scalar ones: A M {x^x 5 ) =^> 
X5), ip(x, x 5 ), which equivalent to infinite number of four-dimensional fields: 

n=oo 

A^ix^xs) = A Kn) a (x)exp(iMnx 5 ) 

n=— 00 
n=oo 

Lp a (x^,x 5 ) = (f(n) a {x)exp(iMnx 5 ) (11) 

n=— 00 

where M ~ A*r n \ = and similarly for (p. Performing integration on x$, one obtains 

four- dimensional gauge theory with infinite dimensional gauge algebra: 

[j. a ± b~i rabcA. c 

\J"n ) J J ^n+m 

In this, gauge transformations of all fields have the form: 

n-k) 

5<p n a = MnA n a + f abc A k b <p (n - k) c (12) 

From the last line we see, that all the scalar fields, except tpo, have non-homogeneous trans- 
formation laws and play the role of Goldstone fields. As a result, such model describes finite 
number of massless vector A^Q a and scalar (fo a fields and infinite tower of massive n / 
vector fields. Introducing appropriate covariant derivatives for these fields: 

D^ n a = d^ n a - / abc A M(fc) V(n- fc) c - MnA Kn) a 
the total Lagrangian could be written in a compact form: 

n=oo y -y 

£ = X! l~jA^ v ( n fA^ u (_ n ) a + -D^{p n a D^{p__ n a ] 

n=—oo 

More complicated examples, corresponding to partial breaking of initial gauge algebra, could 
be found in [47] . 



1.2 Spin 3/2 

Our second example devoted to the possible electromagnetic interactions for massive spin 
3/2 particles P, [10]. As is well known using two fields — vector-spinor ^ and spinor \ one 
can construct gauge invariant formulation of massive spin 3/2 particle. So we start with the 
sum of free massive spin 3/2 and massless spin 1 Lagrangians: 

Tfl — / 3 — 

— + -m{^)x-mxx (13) 
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where gauge transformations leaving this Lagrangian invariant look like: 

5o* M = d^r] + — 7 M ?7, 5 x = y 2 mr /> M m = <^A ( 14 ) 

Here 77 is a spinor parameter, while A is a scalar. We prefer to work with Majorana spinors, so 
in what follows all spinor objects ^ M , x and r\ will be considered as doublets of (anticommut- 
ing) Majorana spinors. In this, minimal gauge interactions corresponds to the replacement 
of ordinary partial derivatives by covariant ones: 

0„ = 0„ + ?4» q ={-e o)' ? 2 = - e2/ ( 15 ) 

As a result of such replacement, the Lagrangian lost its gauge invariance (just because 
covariant derivatives do not commute): 

5 C = MtfA^w, AT = l -e^A aP (16) 

So we try to restore gauge invariance by adding to the Lagrangian non-minimal terms linear 
in gauge field-strength A^ u : 

mC x = ^[ ai A»» + a 2l5 A»" + a 3 g^{aA) + a 4 (A^a a » + a^A a ")]q*„ + 

+i^,(a 5 A^ + a 6l5 Ani»qX + yX?M)x (17) 
as well as all possible linear terms for ^-transformations for all three fields: 

m5i^ M = iq(a 1 A tlu + ol^A^)^ r] mb x \ = qa 3 (aA)r] 

mbxA^ = a^^qr]) +ia 5 (x'J fl qr]) (18) 

Straightforward calculations show that the gauge invariance could be restored (up to terms 
quadratic in A^) provided all the unknown coefficients are expressed in terms of two pa- 
rameters, say «! and a 3 , 



oli = a±, CX4 = 2a±, = —2a 3 , a\ = — 2ai, 

/2 



2 

d2 = 2a\. 0,3 = a>4 = 0, as = clq = — 2«3, a-j = 2\ —a 3 



while this two parameters satisfy the relation: 

2a x - 4^|a 3 + 1 = (19) 

So we have one paramater family of Lagrangians. Now, if we calculate the commutator of 
two r\ transformations on the vector field A^, we obtain: 

[Si, S 2 ]A^ = -Ue\a x 2 + 2a 3 2 )fe7^i)A^ (20) 
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Due to relation on the parameters ot.\ and as it is impossible to set both parameters simul- 
taneously equal to zero. This means that commutator of two rj transformations always give 
a translation and to proceed further on we must introduce graviton as well. So any model 
with minimal gauge interaction for massive spin 3/2 particle in flat space must be a part 
of some (spontaneously broken) supergravity theory. In this case all coupling constants are 
related with gravitational one k ~ — and we obtain usual for supergravity models relation 
betwen spin 3/2 mass, Plank mass and gauge coupling constant m ~ em p i\ Note in pass that 
one more example of constructive approach, namely electromagnetic interactions for massive 
spin 2 particles, could be found in [T5] . 



2 Free massive spin 2 

Our main starting point is the gauge invariant description of free massive spin 2 particle. As 
we will see part of the results obtained crucially depend on the space-time dimension d and/or 
the presence or absence of cosmological term. So in what follows we will work in space-time 
with arbitrary d > 3 dimension and with non-zero (positive or negative) cosmological term. 
We denote metric of such space as g^ u and appropriate covariant derivatives as D^. Till 
the last section, where g^ v will become dynamical, it is just fixed background metric. Our 
convention for the commutator of two covariant derivatives is: 

[Dp, D v \v a = K{g m 8/ - 5/g va )v p (21) 

where k = — ( d _ 1 2 ^_ 2 ) > A — cosmological term. 

One of the nice features of gauge invariant description is that it allows one to consider 
general case of non-zero mass and cosmological term including all possible massless and par- 
tially massless limits [IB]. To describe the massive spin 2 particle in <i-dimensional constant 
curvature space we will use the following Lagrangian: 

-\{D,A U - D u Arf + ^il(D^) 2 + 2m{hTD lt A v - h{DA)) - 



A»D^ '—± >-{h^\ 



d-2 ' ^ 2 v ^ ' 

2(d — l)mc 0l 2d(d — l)m 2 2 2 



where Cq = J m? + n(d — 2). Here symmetric tensor h^ v is the main gauge field, while 
vector A^ and scalar ip are Goldstone fields necessary for gauge invariance. To simplify the 
calculations we have chosen non-canonical normalization for kinetic terms of these fields. 
Also note that there is an ambiguity in the structure of kinetic terms for the h^ u field. 
Indeed, in flat space (where derivatives commute) the second and third terms in the first line 
differ by total divergence only. But in spaces with non-zero cosmological term they lead to 
slightly different structure of mass-like terms and the choice we made gives the most simple 
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and convenient structure of these terms. It is not hard to check that this Lagrangian is 
invariant under the following local gauge transformations: 

2/7Z 

8 Q A^ = D^X + 5 ip = c A (23) 

Let us note, that these and all subsequent results depend crucially on use of canonical 
description of massless spin 2 particles. As is known [19], there exists alternative possibility 
based on traceless symmetric tensor (or on additional Weyl symmetry). Canonical analy- 
sis shows that right number of physical degrees of freedom in this case is achieved due to 
combination of first and second class constraints. So the relation between alternative and 
canonical formulations is similar to that between, say, non gauge invariant and gauge in- 
variant formulations of massive spin 1 particle. Indeed, additional scalar degree of freedom 
(trace h in this case) promotes second class constraints to the first class ones. As we have 
already note in the Introduction, it is hard to formulate constructive approach for theories 
with second class constraints. Moreover, as it was shown in [49j, any attempt to give mass 
to graviton in such alternative formulation unavoidably leads to the appearance of ghosts. 

Note also, that there exists a disagreement on the definition of mass in (A)dS space. 
Indeed, even for the massless particles gauge invariance requires some "mass-like" terms in 
the Lagrangians to be present and one often tries to interpret these terms as real mass. 
Gauge invariant description of massive particles gives, above all, simple and unambiguous 
definition of massless limit, namely it is a limit where Goldstone fields decouple from the 
main gauge field. 

Let us stress the important difference between massless and massive cases for spin 2 par- 
ticle. In the massless case we have one gauge symmetry with vector parameter £ M only, while 
for non-zero mass we have two gauge symmetries with vector £ M and scalar A parameters and 
only these two symmetries together could guarantee the right number of physical degrees of 
freedom even after switching on an interaction. Thus, models based on the £ M invariance 
only like those based on OSP(4) symmetry [50] or models exploring A transformations only 
[5T] always require additional checks for consistency. One more example is a BRST con- 
struction of [52J, where for massive spin 2 particles authors promote gauge transformations 
with vector parameter only (and introduce vector Goldstone field only). As a result, such 
theory describes additional scalar degrees of freedom and all previous results point that it 
must be ghost. At the same time, BRST construction for pure massive spin 2 particle does 
exist [18J, but requires additional scalar gauge transformations and scalar Goldstone field. 

Note also that if one introduces gauge invariant under the A transformations derivatives: 

V^h al3 = - -j—^Arfap, V ^ = D^if - CqA^ (24) 

then the Lagrangian could be rewritten in more simple form: 

c = ^v a ^"v a v - ^v a rv,/i ra - ^(v/i)"(v/i), + (v^^v^ - 

~VW^ - \{D,A V - D v A,f + ^f^(V^) 2 - 

m 2 +K(d-2) 2 2(rf-l)mc t _ 2d{d-l)m 2 2 
2 {k V - M - d _ 2 V + (d _ 2)2 V (25) 
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As could be easily seen from the formulas given above in the massless limit m — > the 
Lagrangian breaks into two independent parts. One of them (for symmetric tensor h^J) gives 
usual description of massless spin 2 particle in (A)dSd background, while the other one (with 
vector and scalar if fields) gives gauge invariant description of massive spin 1 particle 
(or sum of massless spin 1 and spin particles in a flat case). But in de Sitter space k < 
there exist one more special limit cq — > 0. In this, the scalar field ip completely decouples, 
while pair h^, A^ corresponds to the so called partially massless spin 2 particle. Note that 
one can use £ M gauge transformation in order to set vector A^ = 0. In this, the resulting 
simple Lagrangian for h^ v will still be invariant under the A transformations, provided we 
supplement it with restoring gauge A^ = transformation with £ M = — ^D^X: 

1 2m 
5K„ = (D,D U + D V D„)X + -T^9^\ (26) 

We have explicitly checked that such gauge fixed Lagrangian is indeed invariant under this 
transformation. 



3 Self-interaction 

In this section we consider self-interaction of such massive spin 2 particles in linear approxi- 
mation. As we have already mentioned in the introduction, we will not make any suggestion 
on the structure of gauge algebra which could stand behind this theory nor we will not in- 
sist on any geometrical interpretation. Instead, we will construct the most general Lorentz 
invariant cubic terms for the Lagrangian as well as the most general ansatz for gauge trans- 
formations and require that the Lagrangian will be gauge invariant. As a result this section 
will be the most technical part of the paper, but it is important to make sure that our results 
are completely general. 

Even for the massless particles in (A)dSd gauge invariance require introduction of mass- 
like terms into the Lagrangian and appropriate corrections to gauge transformation laws so 
that the structure of the theory resembles that of massive theory in flat space. Thus, working 
with general case of massive particles in (A)dSd, it is convenient to organize the calculations 
just by the number of derivatives. So we start here with the Lagrangian terms with two 
derivatives and gauge transformations with one derivative and analyze all possible vertexes 
compatible with gauge invariance up to the lower derivatives terms. 

/i/i/i-vertex. In this case the most general ansatz for gauge transformations linear in 
/i-field and containing one derivative looks as follows: 

Sh^ = c 1 CD a h llu + c 2 ((Dh)^ u + (Dh) l/ ^) + c 3 g flu (Dh) a C + 

+ c A g^D a hC + c 5 (D^h ua + D u h m )C + c^D^ + D u hQ + 

+ c 7 {h ail D u C + + Csh(D^ + D^) + c 9 VPO + 

+ cio(VA*& + K a D a Q + cug^h^Datf, + c 12 g^h(DZ) 

But due to gauge invariance of free Lagrangian the structure of this transformations could 
be defined only up to arbitrary field dependent gauge transformations (or in other words up 
to possible redefinitions of £ M parameter) which in this case have the form: 

5h, v = d^D^C) + ^(Vf)l + d 2 [D„(h&) + D u (hQ] 
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Moreover, as in all theories where interacting Lagrangian contains the same number of deriva- 
tives as the free one, there always exists a family of physically equivalent Lagrangians related 
by trivial field redefinitions. In the case at hands, we have a possibility to make such redefi- 
nition with four arbitrary parameters: 

hfj,u =^ + Sih^ a h au + s-ihhpv + s^g^li^ 2 + s^g^h 2 

We use all these freedom to bring the gauge transformations to the following more simple 
form: 

<*V = c 1 CD a h llu + c 2 ((Dh)^ l/ + (Dh) u ^) + c 3 g llu (Dh) a C + 
+ c 4 g^D a hC + c 7 {h a ^D u C + h va D^ a ) + c 9 V(^0 

At last, we require that the algebra of these transformations be closed, i.e. 

Simple calculations immediately give c 2 = c 3 = C4 = C9 = and c-j — c±. In the massless case 
the only remaining parameter C\ corresponds to gravitational coupling constant k ~ l/m p i. 
In what follows, we set c\ = 2. Thus, we obtain a very simple final form for these gauge 
transformations: 

5h, u = 2[CD a V + D,CKu + D U CK»] (27) 

In this, 

Certainly, these transformations look exactly like general coordinate transformations for 
covariant second rank symmetric tensor, but let us stress that in our case these are just gauge 
transformations for spin 2 field living in fixed (A)dSd background. Recall also, that concrete 
form of these transformations (i.e. that of covariant tensor and not that of contravariant 
one or tensor density) is just a matter of our choice. This and all our results that follows 
are always defined up to possible field redefinitions of the type shown above. 

Now we construct the most general Lorentz invariant cubic terms with two derivatives: 

Chhh = h^[ ai D^K a (Dh) a + a 2 D l ,K a D a h + a 3 D l ,h af3 D u h a P + ai D^hD u h + 

+ a b D^h al5 D v h aP + a & D^h(Dh) v + a 7 D a h^(Dh) a + a 8 D a h^D a h + 

+ a 9 D a h^D?h au + a 10 D a h^D a h Pu + ail (Dh)^Dh) v ] + h[a 12 D^ h af3 D^h afS + 

+ auDphDOh + a u (DK)^Dhr + a 15 (Dh)^h + a 1& D^h aP D a h^} (28) 

As in the free case, we face an ambiguity because in a flat space not all these terms are 
independent, namely up to total divergence we have (schematically): 

(ai) = (a 5 ) + (a 9 ) - (an) 
(a 2 ) = (a 6 ) + (014) - (ai 6 ) 

so there exists a family of equivalent Lagrangians with two arbitrary parameters. In a 
constant curvature space different choices lead to slightly different structure of mass-like h 3 
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terms and again we use this freedom to bring these terms to the most simple form. Then 
the requirement that the Lagrangian will be gauge invariant (up to lower derivative terms 
given below) gives: 

a i — a n — 3/2, a 2 = —3/2, 03 = —1, 0,4 = = 1, a 5 = 5/2, a$ = —5/2, 

CL-j — Q>lQ = — 2, Og = 2, CLg = d\2 = 1/2; <3l3 = ^14 = 0-16 = — 1/2 

In this, in the non-flat space background we still have non-compensated variations to be 
taken into account in the lower derivatives orders: 

SaChhh + S^o = K{{U-9 > )h al3 D a h^-2{2d-h)h^D tl h aP + 2{d-3)hD ll h- 

- 2(d-3)h ll a D a h+(3d-8)h fl a (Dh) a -2(d-3)h(Dh) ll }^ (29) 

Till now all calculations are the same for massive as well as massless spin 2 particles. 
Let us make here a comment on the so called gravity reconstruction, i.e. on possibility to 
reconstruct full gravity theory starting from free massless spin 2 particle in flat or constant 
curvature space background (53]. As usual, in all calculations of the type given above, all 
the variations which are total divergence are dropped out and of course none of such terms 
can be "reconstructed" in such iterations. But nothing prevent of to add to the Lagrangian 
any such terms to bring the result in a convenient form. As for the ambiguity which arise in 
this iterative process and which is often related with ambiguity in the definition of energy- 
momentum tensor for gravity or matter fields, we have seen that this ambiguity clearly related 
with trivial field redefinitions and up to this freedom the results are essentially unique. 

ZiAA-vertex. Here we start with the most general ansatz for £ M transformations: 

SA„ = c 1 D lx A u e + c 2 D u A lx e + C3(DA)^ + 
+c A A u D^ u + c 5 A v D v ^ + ceA^DZ) 

As in the previous case there exists a freedom connected with the possible redefinitions of 
field A^ and parameter £ M : 

5A^ = d 1 D ll {Ai) A^ -> A^ + sth^A" + s 2 hA„ 

so without lost of generality we can restrict ourselves to 

5 A, = Cl D^A u C + c 2 D u A^C + c 3 (DA)^ 

Usually, then one considers an interaction of gravity with vector fields (abelian or non- 
abelian) one assumes that gravitational field h^ u is inert under the gauge transformations of 
these fields. But now vector field A^ is just a component of massive spin 2 field, so we have 
to consider all the possible gauge transformations with A parameter also. Most general form 
looks like: 

Sh^ = CiD^A y) \ + c 5 g^(DA)\ + cqA^D u) \ + c 7 g flu A a D a \ 
5A^ = CsD^hX + c 9 (Dh)^\ + c 10 hD^X + c u h fny D u \ 
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Once again using the freedom that exists here: 

Sh^ = d 2 D ip (A v) \) 5A^ = dsD^hX) -> + s 3 A^A u + s 4 g^A a 2 

we can reduce these transformations to the form: 

Sh^ = c 5 g^ u (DA)\ 

5A^ = cg{Dh) ll \ + ciohD (l \ + c 1 ih tw D v \ 

Now we proceed by constructing the most general Lorentz invariant cubic terms with two 
derivatives: 

C hAA = h^[ ai D^A u {DA) + a 2 D li A a D v A a + a 3 D^A a D a A u + a 4 D a A^D a A u ] + 
+h[a 5 D a A^D a A l3 + a(i D a ApDpA a + a 7 (DA)(DA)] + 
+a s D a h^D a A tI A u + a 9 (Dhy D ^A a A a + a l0 (Dhf D a A^A a + 
+a ll {Dhy(DA)A^ 1 + a 12 D»hD^A a A a + a 13 D"h(DA)A^ 

Not all these terms are independent here. Indeed, it is easy to check that up to terms without 
derivatives: 

D a h^ u D tI A u A a = h^D a A v D^A a - h^D^A v {DA) + (DhfD^A 

D»hD v A^A u = h{DA)(DA) - hD a A (3 D a A l3 + D fl h{DA)A (1 
D a h fiV D /J ,A a A v = (DhY(DA)A ll + {DKfD^A* - D a h^D^A v A a 

Using this freedom and requiring that the Lagrangian will be gauge invariant we finally get: 

Chaa = c 2 h^A, a A u a - jhA^ 2 (30) 

while the only non-trivial transformation is: 

5A, = c£ v A vlk (31) 

The result obtained is of course familiar and could seems trivial, but for what follows it 
is very important that linear approximation does not fix the value of c 2 coupling constant. 
In the "normal" massless gravity one expects c 2 = 2 but as we will see later on massive 
theory is possible provided c 2 = 1 so that Goldstone field A^ must have non-canonical £ M 
t r ansf or mat ions . 

hh(p-vertex. Now the most general ansatz for £ M transformations has the form: 

Sh^ = d(D^ u + D u (p^) + c 2 g^{D a ^>i a ) + c 3 <p(D^ v + D u ^) + c^pg^D^) 
5<p = c 5 (Dh)^ + c 6 D^he + c 7 h^D^ + csh(D£) 

and using one more time the freedom to make redefinitions: 

Sh^ = d^D^ip^) + D v (<p£n)], -> + s^h^ + s^g^h, y? ^ <p+ s 3 h^h^ v + s 4 h 2 
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one can leave c 2 , C5 and cq as the only non-zero parameters. Then writing most general cubic 
terms: 

C hhv = V [a 1 D^h aP D^h a(3 + a 2 D^D a h^ + a z {DhY{Dh) ll + a i {DhYD lx h + 
+a 5 D^hD^h] + D^aeD^h^ + a 7 D a h^h aP + a^Dhyh^ + 
+a 9 D u hh flu + a w D^hh + On(D/i)^/i] 

and using the fact that up to lower derivative terms we have: 

(a 2 ) = (a 3 ) - (a 7 ) + (a 8 ) 

one can check that there is no non-trivial solution for such vertex. 

AAy>vertex. In this case one has to consider A transformations only, so we write: 

5A^ = ciLpD^X + c 2 D^\ 5ip = c 3 (DA)X + c A A^D^\ 

Possible field and parameter A redefinitions 

6A li = d 1 D ll ((pX), A^ -> + sxtpAp, tp -> tp + s 2 A li 2 

allows one to leave c 3 as the only non-zero parameter. Then considering the most general 
cubic terms: 

C v aa = ^[a l D^A v D^A u + a 2 D^A u D u A^ + a 3 {DA){DA)} + 
+D"( P [a A D li A v A v + a 5 (DA)A li ] 

and taking into account that 

(pD tl A u D u A fl = <p(DA)(DA) + D"<p(DA)A fl - D»D y A^A v 

we obtain the following simple Lagrangian 

£-ipAA = a^A^J 2 (32) 

which is trivially gauge invariant. 

fupip-vertex.. Here the only possible terms for £ M transformations are: 

S(p = c 3 D^e + cM D 

Moreover, using field redefinition 

</?—></? + s\hip 

we can leave c 3 as the only non-zero parameter. Then the requirement that the Lagrangian 

C h99 = dihrDptpDM + a 2 hD^D^ + a 3 (Dhy D ^tptp + a^hD^ip 

will be gauge invariant gives: 

2c 3 (d-l) c 3 (d-l) 
a, = j-^-. a 2 = a 3 = -a 4 
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In this, the arbitrary parameter 03 is related with one more field redefinition (which does 
not change the structure of gauge transformations) 

hua, -> + s^g^ip 2 

So without lost of generality we can set a 3 = and obtain: 

4 W = < ^^[-c 3 h^D^D^+^h(D<p) 2 } 

5ip = c^Dtf (33) 

Again this result could seems trivial but it is important that coupling constant c 3 is not fixed 
yet. 

</?</?<^-vertex. In this last and simplest case we have only one possible term in the La- 
grangian, no non-trivial gauge transformations and one possible field redefinition: 

C v a = diipD^ipD^ip, ip — > ip + siip 2 

showing that this vertex is trivial one. 

Let us collect together all the pieces obtained so far. Total Lagrangian with two deriva- 
tives: 

£0 = Chhh + c 2 h^A fia A u a - -^hAfJ 2 + aotpA/J 2 + 

+ 2 ^^ L [-c 3 h^D,ipD v ip + C jh(D<p) 2 } (34) 
while the only non-trivial gauge transformations look like: 

8A» = c 2 CA^ 8<p = cz£»D ll tp (35) 

Let us stress once again that parameters c 2 and C3 are not fixed yet. 

Now we proceed to the next order. Namely, we construct the most general cubic terms 
with one derivative for the Lagrangian as well as the most general linear terms without 
derivatives for gauge transformations laws. At this order there are no ambiguities related 
with field redefinitions so all calculations are completely straightforward. The resulting part 
of the Lagrangian has the form: 

d = m[W v {Dh) li A v -h{Dh) li A' i + 2h^D^h ua A a - 
-ZhTDphA,, - ZhTDJi^A* + hD^hA^] + 

+2c b h^A^D u p - c b h(ADip) + b^mA^D^pp (36) 

while additional terms to the gauge transformations are: 

4 6 — d 

5h^ u = m[A^ v + A v £n - j^g^iM) + h^X] 

SA^ = mh^f + 2a mip^ Sip = -c c 3 {A£) — ipX (37) 
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Here b = 2c3 J^~^ , a = ^td-l) ■> c 2 = 1 and c 3 = 1 — ^2^-2) ■ ^ s con ^ easily be seen 
from these formulas, in particular from the expressions for the a and C3 parameters, in 
arbitrary space-time dimension d general solution exists for cq 7^ only But in d — 4 
dimensions there exists another solution with cq = then scalar field ip completely decouples 
which corresponds to the self-interaction for partially massless gravity. This fact is clearly 
connected with the conformal invariance of partially massless theories namely in d — 4 
dimensions [37J. Now we proceed by considering general solution with cq 7^ in arbitrary 
space-time dimension but we will comment on partially massless gravity at the end of this 
section. 

The last but not least part of our calculations deals with cubic terms without derivatives 
in the Lagrangian. The most general form for these terms looks like: 

C = b x h! a 'h va h a ' i + b 2 hh» u h tlu + b 3 h 3 + b^hTh^ + b 5 ph 2 + 

+b 6 cp 2 h + b 7 cp 3 + b 8 h^A^A u + b$hA 2 + b w ipA^ (38) 

There are no any new terms for the gauge transformations here, so the gauge invariance 
must be achieved with the gauge transformations obtained at the previous orders. Indeed, 
it turns out to be possible and gives: 

- (7d-16)« , 5m 2 (3d -7k , m 2 (2d-5)« 
b 1 = m 2 + ± — '— b 2 = b 3 = — + y 



b 4 = -2h 



6 4 2 4 6 

2c m(d-l) (d-4)m 3 (d 2 + 5d - 6)m 2 (3d - 2)(d - 4)m 4 



d-2 2c (d-2)' 2(d-2) 2 4c 2 (d-2) 2 

, 9/ ,9 , s d(d + 2)(d - 4)m 5 
b 7 = 2dc 2 (d - lid + 10) + 



6c 3 (d-2) 3 

2 n , 1 ,\ 1 m2 /, ,\ , 2c (d-l)m (d-4)m 3 

6 8 = m 2 - 2(d - 1)«, bg = — + (d-l)K, b 10 = uv ' + K . - 

2 d—2 c (a — 2) 

The dependence of these coefficients on the space-time dimension d and cosmological term 



k (recall that c$ = ym 2 + (d — 2)k) appears to be rather complicated so let us give explicit 
expression for non-derivative terms in d = 4 (any way this dimension is special for us) 

Co = (m 2 + 2K)[h^h ua h afl -^hh^h^ + h 3 ] 

,,,„■, 3c m j9 15m 2 9 , 3m 3 , 

+3c Q m(ph^huu —<ph 2 H ip 2 h y? 3 + 

2 A c 

+ (m 2 - §K)hTApA v + m K hA 2 - ZcompA 2 (39) 
A few comments are in order. 

• One of the most important results of our investigation is that possibility to switch on 
self-interaction in linear approximation crucially depends on the non-canonical form 
of gauge transformations for vector Goldstone field A^ (recall, that results of linear 
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approximation are model independent up to restriction on the number of derivatives). 
Gauge transformations we obtained: 

sh^ = v^e, + + 2(rv a v + D^ a h av + ^rv) 

5 Ay = m^ + CAa^ + mh^C 

resembles very much the situation with massive spin 1 particles which we discussed in 
section 1: 

X A & Pi to, „^abc A b/-c r a „ /-a 9 abc J>£C 

6 Ay, = dpt - ge Ay £ d<p = - -e ip f 

So by analogy with spin 1 case one can suggest that there are two possible ways to 
proceed beyond linear approximation. At first, if one evade an introduction of any 
additional fields in the model, then this results in highly non-linear theory with higher 
and higher derivatives. Indeed, we have explicitly checked that one can not construct 
quadratic approximation without introduction of such higher derivatives terms. On 
the other hand, one can try to introduce additional fields (analog of Higgs field) in 
order to restrict the number of derivatives. 

• Note that exactly as in the free case all the terms containing vector fields (except 
kinetic ones) could be "hidden" if one introduce A-covariant derivatives: 

2m 6 d At, 

v ,j,n a p = Uyh^p — - — ^Ayg aj3 — - — ^mAyh a p 

= D^-coAp + ^y-Atf (40) 

• One can always use the £ M and A local gauge transformations to choose the gauge where 
Ay = and ip = 0. In this, the main difference of massive and massless theories are 
the cubic completion of famous Fierz-Pauli quadratic mass terms: 

m 2 [h^h ua h°"" - ^hh^hy U + hi 3 } 

It is instructive to compare this result with that of [54] (see also [55]) obtained in a 
very different context. The investigation of [54] shows that allowed cubic terms must 
be combination of h^ v h va h a y — hh^ u hy V and h^ u h UQ h a y — ^hh^ u hy U + \h 3 . It is easy to 
see that our terms correspond to combination with coefficients 5, so at this point our 
results agree with that of [54"] . 

• In the discussions of possible mass terms for gravity it is very often assumed that 
the full theory could be the sum of usual massless gravity plus some invariant mass 
terms, e.g. [17]. If we denote effective metric as §y V and assume that in the lowest 
approximation y/— g ~ \f—g + h and g^ v ~ g^ v — c lhF' v then we easily obtain: 

- -y 4~99 m 9 V \K^ - h^Kp) ~ 2m 2 [F/ lra F - -hh^h^ + -h 3 ] (41) 

so the structure of cubic terms does not contradict (up to a factor 2) such assumption. 
Nevertheless, let us stress that in the massive case the gauge transformations of tensor 
hy U is changed so it is hard (if at all possible) to represent massive theory as a sum of 
two parts which are separately gauge invariant. 
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• In the massless limit massive spin 2 particle in (A)dSd decompose into massless spin 
2 particle and massive spin 1 (or massless spin 1 and spin in flat space). In this, 
two gauge transformations £ M and A are completely independent and A transformation 
is just usual gauge transformation of abelian vector field. But then mass m ^ the 
commutator of £ M and A transformations gives: 

[5(A), 5(0] V = m[D^\) + D u {^\)] 

[5(\),5(0]A, = m%\ (42) 

so if we denote a generator of £ M transformation as P M and that of A one as D, we get 
commutation relation for Weyl group [P M , D] ~ mP^. 

• As we have already mentioned, apart form general solution which exists in arbitrary 
dimension d, in dS& and only in this dimension there exist another solution with cq = 
\/m 2 + 2k = which corresponds to self interaction of partially massless spin 2 particle 
in de Sitter background. In this, scalar field if completely decouples, leaving us with 
h^ v and fields (note the absence of h 3 terms): 

C'int = £hhh + A^Ay — —JlA^ + 

+m[4h^(Dh) l ,A l/ - h{Dh)„A<* + 2h /M/ D fl h va A a - 
-ZhrD,JiA v - 3W u D a h^A a + hD^hA^} (43) 

Using local £ M gauge transformation one can always choose the gauge A^ = 0. As in the 
free case, the resulting Lagrangian will still be invariant under A transformation pro- 
vided we supplement it with the appropriate restoring gauge A^ = £ M transformation. 
Indeed, let us look at the form of these transformations when A^ = 0: 

5 V = D^ u + D v £p + mg^X + 2(CD a h^ + D^CKu + D„CK») + mh^X 
6A„ = D^X + m^ + mh^C (44) 
Then it is easy to see that compensating gauge transformation has to be: 

^---[D^X + h^X] 

We have explicitly checked that the Lagrangian for the h^ v field is indeed invariant 
under such combined A and £ M transformations. In the A^ = gauge the Lagrangian 
looks as the linearization of usual gravity in dS± background but this holds in the 
linear approximation only. We have checked that it is impossible to proceed with 
quadratic approximation without introduction of higher derivative terms and/or some 
other fields. 



4 Interaction with matter 

In this section we investigate possible interactions of massive spin 2 particles with matter 
fields of lower spins. The strategy will be the same as in the previous section — we construct 
the most general cubic terms with no more than two derivatives in the Lagrangians and 
corresponding linear terms with no more than one derivative for gauge transformation laws 
with the only requirement that the Lagrangian be gauge invariant. 
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4.1 Spin 



We start with the simplest case of massive spin particles with the usual free Lagrangian: 

Co = \{d^f - ^fn 2 (45) 

All calculations here are very simple and one can easily check that up to possible field 
redefinitions an interacting Lagrangian has a form: 

A = -%[h^d„nd„n-h(d»n) 2 + -<p(d„n) 2 }- 

Z Z Co 

_£^| A _AV (46) 
4 c 

while gauge transformations for it field look like: 

We see that vector field does not couple at all (because ir field is inert under A trans- 
formations), while coupling of scalar field <p clearly shows an ambiguity between flat space 
and massless limits. Indeed, this coupling depends on m/c = mj ^m? + (d — 2)k and we 
have two different limits. On one hand, we can take a massless limit m — > while keeping 
cosmological term k small but non-zero. In this, scalar field <p completely decouples so that 
massless limit of massive theory agrees with purely massless theory results. On the other 
hand, if we take flat space limit k — > while keeping non-zero mass m, then the coupling 
constant for scalar field (p tends to 1 and does not depend on mass any more. As a result, 
in the massless limit scalar field ip does not decouple from matter field ir. 



4.2 Spin 1 

Our next example — interaction with massive spin 1 particles. Due to our usage of gauge 
invariant description this includes massless limit as well. Thus, we introduce two fields - 
vector and Goldstone scalar n and start with the free Lagrangian which has its own gauge 
invariance: 

Co = -\b, u 2 + \{D^f - mi B»D^ + ^B, 2 
SB^ = c^A, 57i = m 1 \ (47) 

An investigation of possible hBB and Innr vertexes with two derivatives goes exactly the 
same way as that of hAA and htptp vertexes in the previous section, so we will not repeat 
these calculations here and write the corresponding part of Lagrangian: 

A = C ±WB m B va - jhB, v 2 - ^hTD^Dv* + jHD^f + 

+a ip(D^) 2 + a^B^ 2 + a^A^B^ (48) 
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which is invariant under the following transformations: 

5B„ = Cl CB u „ 5n = c 2 ^D^ (49) 

Now we proceed by adding all possible lower derivatives terms to the Lagrangian and corre- 
sponding terms to gauge transformations. This procedure results in: 

1 in 

C 2 = m lCl [h^B^D u Ti - -h(BD-n) + — V (BDtc)] - 

2 Co 

_^2L [hrBitBv _ l hB * + ™^ 2] (50) 

2 2 Co 

with the only new term in gauge transformations: 

S'tt = -dm^B^ 

In this, gauge invariance requires c\ = c 2 , a = ^, a\ = ^^Z^) ' a 2 = 0. Once again, our 
results show an ambiguity betwen flat space and massless limits. Only for massless vector 
field in d — 4 dimensions this ambiguity is absent. 



4.3 Spin 1/2 

Our last example — interaction with massive spin 1/2 particles (recall that we prefer to work 
with Majorana spinors). As is well known, to describe spinor fields living in curved back- 
ground one has to use first order formalism in terms of "tetrad" e^ a and Lorentz connection 
so that 7 M = e /J(l 7 a and 

1 K 
[D„, D U ] X = -R^ ab cr a bX = 2°>* ( 51 ) 

In this, to describe an interaction of spinor fields with our massive spin 2 particles it is also 
convenient to use first order formulation of such particles [56J in terms of three pairs of fields: 
(e Ma , oj fl ab ), (A^,F ab ) and (ip,7i a ). We will not reproduce Lagrangian for such formulation 
here (it could easily be found in [56]) because all we need here is the structure of gauge 
t r ansf or mat ions : 

r j I 

5h^ a = D^ a + Tj^ a + — e^ a X 

5u, ab = D^-^ie^-eX) (52) 
<L4 M = D^X + 5ip = c A 

where apart from £ a and A gauge transformations we have now one more transformation 
with parameter r} a b = —r]b a . Again we start with the free Lagrangian for massive spinor \: 

C = -x^D^x - ^XX (53) 

and construct the most general cubic terms (this time with no more than one derivative) 
compatible with all gauge symmetries. We obtain for interaction Lagrangian: 

_»*/»(* ™ )a (54 ) 



2 v 2c 
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while for gauge transformations we get: 

5X = eD, X - -^r))x - —Ax (55) 

In this case an ambiguity between flat space and massless limits exists for massive mi/2 7^ 
spinor field only, while the results for massless spinor agree with that of purely massless 
gravity. 

5 Interaction with gravity 

There exist well known difficulties one faces in any attempt to construct interacting theory 
for a collection of massless spin 2 particles [38l [391 EH SD H2]. For example, for the case of 
just two massless spin 2 particles there are only two possibilities. One of them corresponds to 
two copies of usual gravities completely independent of each other, while the other possibility 
which does have interaction requires that one of these particles has wrong sign of kinetic 
terms and be a ghost. Moreover, there are examples of theories of such kind, where one of 
spin 2 particles is massive, coming from higher derivative gravity models e.g. [57]. Also, as 
is well known, there are examples of consistent theories with infinite number of massive spin 
2 particles in Kaluza-Klein models, but it could be shown [15J HH] that it does not contradict 
with general results of [38l [39] . 

But all this still leaves a possibility to construct a consistent theory where one massless 
spin 2 particle interacts with one or several massive spin 2 ones. In this section we investigate 
the case of one massless and one massive particles. In this, there are two possible cubic ver- 
texes (besides self-interaction). We have checked (though we will not present details of these 
calculations here) that it is impossible to construct linear approximation with two massless 
and one massive particles (once again up to the restriction on the number of derivatives). 
As for the other case, i.e. interaction of massive spin 2 particles with usual gravity, it is not 
hard to constract general covariant vertex with arbitrary number of massless fields but still 
bi-linear in massive fields, as we are going to demonstrate. 

First of all, let us note that in this section metric g^ is not just a fixed background any 
more, but it is a dynamical field with its own equation of motion: 

1 K(d - l)(d - 2) , . 

R^u - -g,uR + f 9yv = (56) 

As usual in gravity, we will assume that connection is metric compatible D a g^ u = and we 
have usual identities: 

D»R^ = X -D V R 

As in the case of massive particle living in constant curvature background, it is convenient 
to organize the calculations by the number of derivatives. So we start with the sum of kinetic 



20 



terms of our three fields h^ u , and ip: 

C 2 = ^D a h^D a h^ v - D a h^D^h ua + {DhyDfJi - l -D»hD^h - 

~\{D,A V - D v A,f + ^L^lFipDpip (57) 

and corresponding gauge transformations with one derivative: 

Si V = D £* + D ^ M„ = D ^ ( 58 ) 

Covariant derivatives do not commute and, as a result, this Lagrangian is not invariant under 
these gauge transformations. But gauge invariance could be restored (up to lower derivative 
terms we will take into account later) if one adds to the Lagrangian: 

AC = -2R^h^ a h ua + BTh^h + ^Rh afS h a p - ^Rh 2 (59) 

and requires that metric g^ u has non-trivial transformation 

5 l9aP = 2{DJi aP - D a hfr - D p h aix )e (60) 

Now we proceed with lower derivative terms. The part of the Lagrangian with one derivatives 
turns out to be: 

A = 2m[h" v D^A v - h(DA)\ - ^~^ A»D^ (61) 

while corrections to gauge transformation look like: 

2m 

8oh» v = _ ^ g^X, 5 A^ = m^, 5 (f = c A (62) 

Also this requires additional modification of g^ u gauge transformations: 

2 2m(d — 4) 
^9 a p = 2m[A a ^ + Apta - _ 2 ^ 9ap(M)] [d-2) hapX ^ 

At last, we must add possible terms without derivatives, the most general form being: 

C = —h^h^ + |/i 2 + b 3 h<p + b^ 2 + b b A 2 (64) 

There are no any additional corrections to gauge transformations at this level, so the gauge 
invariance must be achieved with the ones we already have. And indeed it turns out to be 
possible giving: 

b 1 = m 2 + K (d -l)(d- 2), b 2 = m 2 + -(d - l)(d - 2) 
2(d-l)mc 2m 2 d(d-l) 

h = — (J^T' 64 = ( d -2y ' 65 = Md ~ 1} 
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This time also there is a possibility greatly simplify all the expressions by introduction of 
A-covariant derivatives: 



2m 




(65) 



(d-2) 



In this, total Lagrangian could be rewritten in the form: 



C 




(66) 



while gauge transformations for the metric field g^ u look like: 



5g a/3 = 2(V M /i a/3 - - Vphan)^ + 



2m(d - 4) 
(d-2) 



(67) 



Thus, in this approximation an interaction of massive spin 2 particles with gravity exists 
in any space-time dimension d > 3 (though from the last equation one can see that d = 4 
case is also special here). In particular, nothing prevent us to consider the c — > limit, 
i.e. interaction of partially massless spin 2 particle with gravity. But as in the case of 
self-interaction, as we have explicitly checked, if one tries to proceed with terms quartic in 
massive fields than one will find that higher derivatives interactions and/or some additional 
fields are necessary. It is instructive to compare our results here with the investigations 
of massive spin 2 particle in gravitational background [TTJ [121 H3]- I n general, results are 
similar, but the structure of Rhh terms is slightly differrent. 

Conclusion 

We hope that the main lesson from our paper is that constructive approach based on gauge 
invariant description of massive high spin particles does allow one efficiently investigate 
possible interactions of such particles. It is important that due to peculiarity of linear 
approximation the results obtained for any particle are completely model independent and do 
not depend on the presence of any other fields in the theory. In particular, an impossibility 
to construct an interaction in linear approximation means that such interaction does not 
exist at all. One of the examples of such "no-go" results is the absence of self-interaction for 
partially massless spin 2 particles in d ^ 4 dimensions. 

One of the striking and unexpected results is that the existence of self-interaction for 
massive spin 2 particles crucially depends on non-canonical form of gauge transformations 
for Goldstone A^ field. By analogy with massive spin 1 case, one can assume that in full 
interacting theory (if it exists at all) one must deal with non-linear realization of £ M symmetry 
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with higher and higher derivatives. We are not aware on any works on non-linear realization 
of space-time symmetries where something similar appears. 

In the investigation of the massive spin 2 particle interacting with matter (i.e. spin 0, 1, 
and 1/2 particles) out results clearly show the ambiguity betwen flat space and massless limits 
which reveals itself through the dependence of scalar Goldstone field (p coupling constant on 
mass and cosmo logical constant. Let us stress once again that this results are also model 
independent. 

Throughout the paper we restrict ourselves with interaction terms with no more than two 
derivatives in the Lagrangians (and correspondingly no more than one derivative in gauge 
transformations). But many of our results clearly show that to construct full interacting the- 
ory beyond linear approximations one unavoidably will have to introduce higher derivatives 
interactions. But such higher derivatives interactions could, in principle, change the results 
obtained here for linear approximation. This question deserves further study. 
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